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1.  INTRODUCTION 


The  work  presented  here  is  a  continuation  of  attempts  to  unravel  details  of  reactions  that  occur  in  the 
decomposition  of  hexahydro-l,3,5-trinitro-l,3,5-triazine,  more  commonly  known  as  RDX.  The  primary 
decomposition  of  this  nitramine  under  collisionless,  infrared  multiphoton  excitation  (IRMPE)  conditions 
occurs  through  two  competing  reactions  (Zhao,  Hintsa,  and  Lee  1988),  denoted  (PI)  and  (P2)  below.  The 
reactions  are  scission  of  the  weakest  bond  (N-N02)  to  form  an  unstable  cyclic  radical  and  N02: 

RDX  — >  C2H4N404  +  N02  (PI) 

and  concerted  decomposition  of  the  ring  to  form  three  highly  vibrationally  excited  methylene  nitramine 
(MN)  molecules: 

RDX  ->  3CH2NN02.  (P2) 

These  reactions  were  determined  by  analyzing  time-of-flight  (TOF)  mass  spectra  of  products,  and 
comparing  these  with  predictions  based  on  several  proposed  decomposition  schemes  (Zhao,  Hintsa,  and 
Lee  1988).  Reactions  that  could  not  explain  various  components  of  the  TOF  spectra  were  then  eliminated 
from  the  series  of  proposed  reactions.  In  the  case  of  the  primary  decomposition  of  RDX,  reactions  (PI) 
and  (P2)  were  the  two  reactions  that  were  consistent  with  the  experimental  measurements,  particularly 
concerning  the  TOF  spectra  of  the  large  mass  products  (m/e  >  74).  However,  a  series  of  secondary 
decomposition  reactions  of  RDX  had  to  be  invoked  to  fit  all  of  the  observed  TOF  spectra  (Zhao,  Hintsa, 
and  Lee  1988). 

The  secondary  reactions  of  (P2)  proposed  by  Zhao,  Hintsa,  and  Lee  (1988)  are  all  concerted  molecular 
eliminations.  They  did  not  observe  evidence  of  the  simple  scission  reaction  of  the  weakest  bond  in  MN: 

CH2NN02  H2CN  +  N02.  (I) 

They  propose  two  channels  for  the  decomposition  of  MN  (Zhao,  Hintsa,  and  Lee  1988): 

CH2NN02  ->  HCN  +  HONO  (or  HN02)  (II) 

and 
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(Ill) 


ch2nno2  -»  n2o  +  h2co. 

The  TOF  spectra  indicate  that  OH)  is  a  minor  channel  compared  to  the  decomposition  of  MN  through  (II). 
The  results  of  the  experiments  did  not  determine  whether  HONO  or  HN02  is  the  more  probable  isomer 
of  (II). 

There  have  been  suggestions  that  the  cyclic  radical  C>H4N404  of  (PI)  spontaneously  decomposes  to 
form  MN  (Schroeder  1985a;  1985b).  Regardless  of  the  parentage  of  MN,  the  determination  of  its  fate  is 
important  in  understanding  RDX  decomposition.  The  purpose  of  this  study  is  to  examine,  by  using 
classical  trajectories,  the  reactions  of  highly-excited  CH2NN02,  which  may  be  a  product  of  RDX 
decomposition.  Unfortunately,  the  experimental  measurements  do  not  provide  detailed  information  about 
the  internal  energy  distributions  of  the  nascent  CH2NN02.  For  example,  it  is  not  known  whether  the 
energy  is  distributed  statistically.  Zhao,  Hintsa,  and  Lee  (1988)  estimate  that  the  average  internal  energy 
of  the  RDX  following  IRMPE  is  approximately  80  kcal/mol.  This  information,  in  conjunction  with  a 
Hess’  law  analysis  of  the  endothermicity  of  (P2)  provided  by  Sewell  and  Thompson  (1991),  suggests  that 
the  MN  molecules  formed  from  (P2)  have  total  energies  ranging  from  55-65  kcal/mol.  This  is  sufficient 
energy  for  MN  to  undergo  the  decomposition  reactions  (I)— (III).  Classical  simulations  of  the  primary 
decomposition  of  RDX  (Sewell  and  Thompson  1991)  show  that  most  of  the  energy  in  the  nascent 
CH2NN02  is  vibrational. 

Methylene  nitramine  has  never  been  isolated.  The  only  characterizations  of  this  molecule  come  from 
theoretical  predictions.  In  a  previous  report,  we  developed  a  potential  energy  function  of  this  molecule 
(Rice  et  al.,  submitted  for  publication).  We  based  the  potential  energy  surface  (PES)  of  MN  that  we  used 
in  this  work  on  the  ab  initio  multiconfigurational  (MC)  SCF  and  multireference  (MR)  Cl  electronic 
structure  calculations  of  Mowrey  et  al.  (1990;  unpublished).  They  calculated  the  structure,  relative 
energies,  and  frequencies  of  equilibrium  MN,  the  five-centered  transition  state  for  HONO  elimination,  and 
the  asymptotes  for  the  N-N  scission  reaction  (Mowrey  et  al.  1990).  Additionally,  they  calculated  points 
along  the  reaction  path  from  the  five-centered  transition  state  of  (II)  leading  to  reactants  (MN)  and 
products  (HONO  +  HCN)  (Mowrey  et  al.  unpublished).  After  adjusting  for  zero-point  energy  effects, 
basis  sets  and  different  levels  of  theory,  they  estimate  that  the  activation  energy  for  (II)  is  31  ±4  kcal/mol. 
They  also  estimate  the  N-N  bond  dissociation  energy  for  (I)  to  be  35  ±  4  kcal/mol.  Although  these 
estimates  suggest  that  reaction  (II)  is  energetically  favored,  the  steric  effects  associated  with  reaction  (II) 
would  make  reaction  (I)  entropically  favored. 
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Mowrey  et  al.  (1990)  suggested  that  because  Zhao,  Hintsa,  and  Lee  (1988)  did  not  observe  evidence 
of  reaction  (I),  either  (a)  the  microcanonical  rate  constant  for  (I)  has  an  anomalous  energy  dependence 
which  makes  it  not  probable,  or  (b)  the  barrier  height  for  (II)  is  substantially  lower  than  (I)  and  the  nascent 
MN  do  not  have  enough  energy  for  (I)  to  occur.  Mowrey  et  al.  (1990)  have  investigated  the  latter  and 
have  shown  that  the  activation  energy  for  (II)  is  lower  than,  or  at  least  comparable  with,  the  N-N  bond 
dissociation  energy.  In  the  present  study,  we  have  used  detailed  dynamical  modeling  to  investigate 
whether  MN  exhibits  an  anomalous  energy  dependence  of  the  decomposition  rates.  We  have  calculated 
the  microcanonical  rate  constants  for  reactions  (I)  and  (II)  as  functions  of  energy  by  using  classical 
trajectories  computed  on  an  analytical  PES  that  is  based  on  the  Mowrey  et  al.  ab  initio  results  (1990; 
unpublished).  Additionally,  mechanistic  details  of  the  reactions  are  presented. 

2.  POTENTIAL  ENERGY  SURFACE 

The  analytic  PES  that  we  used  in  this  study  is  described  in  (Rice  et  al.,  submitted  for  publication). 
Structures  corresponding  to  the  reactant  and  transition  state  agree  with  the  ab  initio  values  (Mowrey  et  al. 
1990)  to  within  0.5%.  The  normal  modes  of  vibration  for  MN  are  in  good  agreement  with  the  ab  initio 
calculations  (Mowrey  et  al.  1990).  Also,  the  analytic  function  (Rice  et  al.,  submitted  for  publication)  was 
fitted  to  points  along  the  reaction  path  leading  to  products  for  (II)  and  are  in  good  agreement  with  the  ab 
initio  values  (Mowrey  et  al.  1990).  The  model  PES  predicts  an  activation  energy  for  (II)  of  31.8  kcal/mol, 
and  the  dissociation  energy  of  the  N-N  bond  is  35  kcal/mol  (Rice  et  al.,  submitted  for  publication),  in 
agreement  with  the  Mowrey  et  al.  estimates  (1990). 

3.  METHODS 

Trajectories  at  six  energies  were  integrated  using  a  variable-step  Adams-Moulton  predictor-corrector 
integrator  (Miller  and  George  1972),  with  relative  error  tolerance  set  to  10-7.  Initial  conditions  were 
selected  as  follows:  The  molecule  was  initially  in  the  equilibrium  geometry.  The  total  energy  of  the 
system,  including  the  zero-point  energy,  was  equipartitioned  among  the  normal  modes  in  the  form  of 
kinetic  energy.  A  warmup  trajectory  of  approximately  0.2  ps  was  integrated,  allowing  a  partial 
repartitioning  between  potential  and  kinetic  energy.  These  points  in  phase  space  were  then  used  as  starting 
points  for  Metropolis  Monte  Carlo  sampling  (Metropolis  et  al.  1953;  Raff  and  Thompson  1985).  An  initial 
markov  walk  of  106  steps  was  taken  to  ensure  that  the  system  was  relaxed  away  from  the  initial 
configuration.  A  trajectory  was  then  integrated  until  the  end  tests  indicated  a  reaction  had  occurred,  or 


3 


until  the  simulation  time  exceeded  20  ps.  A  sequence  of  4,000  markov  moves  was  taken  from  the  starting 
point  of  the  previous  trajectory,  and  the  integration/m arkov  walk  pattern  was  repeated  until  an  ensemble 
of  at  least  2,000  trajectories  was  generated.  Three  heavy  atoms  and  one  hydrogen  atom  were  moved  for 
each  markov  step.  Step  sizes  for  the  markov  walk  (given  in  Table  1)  were  selected  to  achieve  50% 
acceptance  of  the  attempted  moves.  Additionally,  the  maximum  bond  extensions  were  limited  in  the 
Metropolis  sampling  to  confine  the  sampling  to  reactant  configuration  space.  The  N-N  and  C-H  bonds 
were  limited  to  less  than  2.0  and  1.2  A,  respectively,  and  the  H-0  bond  is  restricted  to  be  greater  than 

1.6  A. 


Table  1.  Parameters  Used  in  Metropolis  Sampling 


Atom  Type 

Maximum  Coordinate 
Displacement  (A) 

Maximum  Momentum 
Change  (amu-A/t.u.a) 

H 

0.035 

0.025 

C 

0.028 

0.025 

N 

0.025 

0.028 

O 

0.025 

0.035 

a  1  Lu.  =  1.018066  x  10"14  s. 


The  end  tests  used  to  identify  the  reaction  type  are  listed  in  Table  2.  There  is  sufficient  energy  in 
these  simulations  to  allow  for  secondary  reactions  following  the  primary  reactions  (I)  and  (II).  The 
secondary  reactions  observed  in  the  simulations  result  from  decomposition  or  isomerizations  of  the  nascent 
HONO: 


HONO'  — »  HO  +  NO' 

(NO  bond  scission) 

(HI) 

HONO'  ->  ONO'H 

(Hydrogen  Atom  Migration) 

(IV) 

HONO'  ->  ONO'  +  H 

(Hydrogen  Atom  Elimination) 

(V) 

For  purposes  of  discussion,  the  two  identical  oxygen  atoms  listed  in  (III)-(V)  are  distinguished  as  O  and 
O'.  We  will  use  these  same  distinctions  later  in  the  discussion  of  the  results.  If  a  trajectory  triggers  the 
end  tests  specifying  HONO  +  HCN  formation,  ends  tests  for  the  secondary  reactions  are  invoked.  At  this 
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Table  2.  End  Tests 


Primary  Decomposition  Reactions: 

ch2nno2  h2cn  +  no2 
•  Rjvjjvj  >  5.0  A,  Rcn(i)  ^  rch(2)  ^  2.0  A 

CH2NN02  ->  HONO'  +  'HCN 

•  RcH  >  2.0  A,  Rqh  <  1.393  A,  Rch-  <  2.0  A  and  D^,  DCH  <  10-4  eV 

Secondary  Decomt>osition/Isomerization  Reactions: 

HONO'  -»  H  +  ONO' 

•  Rj^q  >  5.0  A 

HONO'  ->  ONO'H 

*  RHO  >  RHO' 

HONO'  -»  HO  +  NO' 

*  rno  ^  rno'  ^  2-0  A 


point  in  the  trajectory  calculation,  the  integration  continues  up  to  the  maximum  integration  time  (20  ps), 
or  until  the  O-H  or  N-0  bond  distances  exceed  5.0  A. 

The  lifetime  of  a  trajectory  resulting  in  (I)  was  taken  to  be  the  integration  time  up  to  the  last  inner 
turning  point  of  the  N-N  bond  vibration.  The  lifetime  of  a  trajectory  resulting  in  (II)  was  taken  to  be  the 
integration  time  up  to  the  point  at  which  all  end  tests  for  reaction  (II)  were  satisfied. 

4.  RESULTS  AND  DISCUSSION 

Table  3  is  a  summary  of  trajectory  results  for  total  energies  (including  the  zero-point  energy, 
27.5  kcal/mol)  ranging  from  62  to  131  kcal/mol.  At  least  2,000  trajectories,  including  those  recounted 
according  to  standard  Metropolis  sampling  procedures  (Metropolis  et  al.  1953;  Raff  and  Thompson  1985), 
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Table  3.  Ensemble  Results 


Primary  Decomposition  Reactions 

Energy8 

(kcal/mol) 

Total  Trajectories 

Total  Reactions15 

Reaction  Ic 

Reaction  IIC 

62.1 

4635 

364  (8) 

322  (88) 

42  (12) 

73.6 

1997 

787  (39) 

673  (86) 

114  (14) 

85.1 

1995 

1646  (82) 

1274  (77) 

372  (23) 

96.7 

2001 

1968  (98) 

1523  (77) 

445  (23) 

|  108.2 

1963 

1963  (100) 

1340  (68) 

|  119.7 

1946 

1946  (100) 

1199  (62) 

747  (38) 

131.3 

1938 

1938  (100) 

1145  (59) 

793  (41) 

Secondary  Decomposition  Reactions  of  HONO  || 

Energy8 

Unreactedd 

HO+NO+HCNd 

H+ONO+HCN^ 

H-Atom  Migration 

(kcal/mol) 

HONO+HCN 

in  HONO 

62.1 

42  (100) 

0 

4 

73.6 

114  (100) 

0 

2 

85.1 

363  (98) 

9(2) 

0 

10 

96.7 

421  (94) 

24  (6) 

0 

10 

108.2 

585  (94) 

38  (6) 

0 

27 

119.7 

663  (89) 

84  (11) 

0 

84 

131.3 

654  (82) 

130(16) 

9(1) 

118 

8  Includes  the  zero-point  energy. 

Values  in  parentheses  are  percentage  of  trajectories  resulting  in  reaction. 
^  Values  in  parentheses  are  percentage  of  the  total  reactions. 

Values  in  parentheses  are  percentages  of  reaction  II  totals. 


were  used  in  ensemble  averaging  for  each  energy.  Because  of  the  low  number  of  reactions  at 
62.1  kcal/mol,  approximately  4,600  trajectories  were  calculated  to  attain  better  statistics.  Trajectories  that 
satisfied  the  end  tests  for  N-N  bond  scission,  but  had  not  undergone  at  least  one  N-N  vibration  were  not 
included  in  the  ensemble  averages  or  reaction  totals.  The  number  of  these  rejected  trajectories  increased 
with  increasing  energy,  leading  to  the  fewer  trajectories  available  for  averaging  at  the  higher  energies. 

The  numbers  in  parentheses  in  the  column  entitled  "Total  Reactions"  are  the  percentages  of  CH2NN02 
that  react  within  the  integration  time  of  20  ps.  At  energies  greater  than  97  kcal/mol,  every  trajectory 
results  in  reaction  within  this  time. 
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Table  3  also  lists  the  numbers  of  trajectories  (and  reaction  percentages  in  parentheses)  resulting  in 
reactions  (I)  or  (II)-  The  majority  of  reacting  trajectories  over  the  energies  sampled  result  in  reaction  (I), 
but  this  majority  decreases  at  higher  energies  as  reaction  (II)  becomes  more  probable.  Additionally,  there 
is  an  increase  in  the  secondary  N-0  scission  and  H-atom  migration  reactions  of  HONO  as  energy 
increases.  The  increase  in  both  of  these  secondary  reactions  as  the  total  energy  increases  is  not  surprising. 
The  increase  in  the  secondary  reactions  of  HONO  is  probably  due  not  only  to  the  details  of  HONO 
formation,  but  also  pn  the  internal  energy  distribution  of  the  CH2NN02  before  reaction.  As  the  total 
formal  energy  of  CH2NN02  increases,  the  energy  distributed  among  the  internal  modes  increases, 
including  the  vibrational  modes  corresponding  to  the  ONO  moiety.  Upon  reaction,  additional  energy  is 
released  to  the  fragments.  Because  the  ONO  moiety  does  not  undergo  bond  breaking  or  significant 
geometric  change  in  the  transformation  from  CH2NN02  to  HONO,  it  is  reasonable  to  assume  that  the 
ONO  moiety  maintains  (or  even  gains)  vibrational  excitation,  some  of  which  result  in  the  secondary 
reactions  corresponding  to  either  N-0  bond  scission,  or  the  large-amplitude  bending  motion  of  HONO 
to  accomplish  H-atom  migration. 

We  did  not  observe  enough  secondary  HONO  decomposition  reactions  resulting  in  H  +  ONO  to 
deteimine  the  energy  dependence  of  this  channel.  We  observed  these  reactions  only  at  one  energy,  and 
the  mechanism  appears  to  be  the  same  for  each  event.  The  mechanism  is  illustrated  in  Figure  1.  Figure  1 
shows  the  behavior  of  the  N-0,  N=0',  and  H-0  distances  during  one  of  the  trajectories  that  results  in  H- 
atom  elimination  from  HONO.  In  this  trajectory,  HONO'  +  HCN  forms  within  0.1  ps.  Upon  formation 
of  HONO,  the  nascent  H-0  bond  is  highly  vibrationally  excited.  The  energy  localized  in  this  bond 
transfers  to  the  remaining  vibrational  modes,  including  the  adjacent  N-0  bond.  The  N-0  bond  appears 
to  experience  increasing  excitation  up  to  3.0  ps,  evident  by  the  large-amplitude  motion.  At  this  time,  the 
energy  appears  to  flow  out  of  the  N-0  bond  and  back  into  the  adjacent  H-0  bond,  also  evident  by  the 
large-amplitude  oscillations  of  this  bond  at  times  greater  than  3.0  ps.  At  3.75  ps,  enough  energy  has 
transferred  back  into  the  H— 0  bond  to  cause  dissociation  to  H  +  ONO .  Throughout  the  trajectory,  the 
terminal  N=0'  bond  does  not  experience  significant  vibrational  excitation. 

A  more  typical  trajectory  resulting  in  secondary  decomposition  of  HONO  to  diatomic  products  is 
illustrated  in  Figure  2.  This  trajectory  also  shows  a  hydrogen  migration  reaction  in  the  HONO  before  the 
decomposition  event.  At  approximately  0.6  ps,  CH2NN02  form  HO'NO  +  HCN,  which  is  most  evident 
by  the  behavior  of  the  H-0'  bond  in  Figure  2d.  After  approximately  1  ps,  the  hydrogen  migrates  on  the 
HO'NO  molecule  to  form  HONO'.  After  another  4  ps,  the  N-0  bond  breaks  and  the  diatomics  HO  +  NO' 
are  formed. 


7 


3.0  -i 


0 


i 


12  3 


TIME  (ps) 

Figure  1.  Internal  coordinates  (a)  N-0  and  N-0'  and  (b)  H-0  in  HONO'  formed  from  primary 
decomposition  of  CH^NNOo  during  a  trajectory  that  results  in  secondary  H-0  bond  scission 
of  the  HONO'  molecule. 
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Figure  2.  Internal  coordinates  (a)  N-O.  fhl  N-Q'.  (c)  H-O,  and  (d)  H-Q/  in  HO'NO  formed  from  primary 
decomposition  of  CFUNNO.  during  a  trajectory  that  results  in  hydrogen  migration  followed  by 
N-Q'  scission  of  the  HO'NO  molecule. 
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Table  4  lists  the  decay  coefficients  and  branching  ratios  for  the  primary  decomposition  reactions  of 
CH2NN02  as  functions  of  energy.  For  a  system  of  molecules  undergoing  unimolecular  decomposition 
through  two  competing  channels,  the  overall  rate  of  decomposition  can  be  expressed  as: 

kTotal  =  ki  +  kn  (!) 

and  can  be  determined  from  trajectory  results  by  fitting  lifetimes  to: 

ln[P]  =  -kTotalt ,  (2) 

where  [P]  is  the  fraction  of  unreacted  CH2NN02  at  time  t.  The  individual  rate  coefficients,  kx  and  kn, 
can  be  extracted  from  the  branching  ratio  (defined  in  this  calculation  as  the  ratio  Nj/Nn,  the  total  number 
of  N-N  bond  scission  reactions  to  the  total  number  of  HONO  formations)  using  the  following  relation: 

Ni/Nn  =  kj/kjj  .  (3) 

The  decay  curves  and  corresponding  branching  ratios  are  shown  in  Figure  3.  All  of  the  decay  curves  and 
branching  ratios  are  linear,  indicating  that  the  rate  for  this  system  is  first-order  and  time  independent. 
Linear  least-squares  fits  of  Equations  (2)  and  (3)  to  these  curves  resulted  in  the  rate  coefficients  listed  in 
Table  4. 

We  have  also  fitted  the  total  decay  rate  and  individual  rates  as  a  function  of  energy  to  the  RRK 
statistical  model  (Robinson  and  Holbrook  1972) 


k(E)  =  v 


"(s-l) 

t 

) 


(4) 


where  E0  is  set  to  35  kcal/mol. 

The  results  of  the  fit  to  the  three  rate  coefficients,  kTo[a],  kx,  and  kn,  as  a  function  of  energy  are  shown 
in  Figures  4a  and  4b.  Values  of  the  best-fit  parameters  are  given  in  Table  4.  Because  this  system  is  so 
well-behaved  and  kTotal,  kx,  and  kn  are  well-described  by  the  RRK  model,  we  were  able  to  obtain  a 
functional  description  of  the  branching  ratios 
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The  branching  ratios  and  prediction  using  the  RRK  parameters  provided  in  Table  4  are  shown  in 
Figure  4c.  If  the  behavior  of  this  system  can  be  well-described  by  this  model  for  energies  much  greater 
than  those  sampled,  then  the  limit  at  high  energy  of  the  branching  ratio  is  0.47,  indicating  that  most  of 
the  high-energy  reactions  are  concerted.  The  results  indicate  that  the  microcanonical  reaction  rates  have 
no  anomalous  energy  dependence  and  are  very  well-behaved.  This  system  does  not  exhibit  nonstatistical 
behavior  for  microcanonical  initial  conditions,  which  ensures  a  statistical  distribution  of  total  energy 
among  the  internal  modes  of  the  molecule. 


Table  4.  Decay  Coefficients  and  Branching  Ratios 


Energya 

kTotal  (PS_1) 

kj  (ps-1) 

kn  (ps  *) 

kj/kn 

62.1 

0.004 

0.0036 

0.00037 

9.6 

73.6 

0.021 

5.7 

85.1 

0.070 

0.019 

3.7 

96.7 

wmmmm 

0.18 

0.053 

3.4 

|  108.2 

0.40 

0.27 

0.13 

2.1 

119.7 

0.68 

0.42 

1.6 

131.3 

1.07 

0.61 

HEEHI 

1.3 

v  (ps"1) 

42.9326 

13.7587 

29.1739 

— 

12.915 

10.9636 

14.6706 

— 

Table  5  lists  the  average  product  translational  and  rotational  energies  for  the  various  primary  and 
secondary  decomposition  reactions  of  CH2NNO2.  The  numbers  in  parentheses  are  the  percentages  of 
energy  available  to  the  products  in  translation  and  rotation.  For  reaction  (I),  most  (-60%)  of  the  available 
energy  is  in  vibration,  with  slightly  more  energy  going  into  relative  translation  than  into  the  rotational 
modes.  The  results  listed  for  reaction  (II)  are  compiled  from  those  trajectories  that  did  not  undergo 
secondary  decomposition  of  HONO.  In  contrast  to  reaction  (I),  approximately  1/3  of  the  available  energy 
goes  into  relative  translation.  HCN  experiences  almost  double  the  rotational  excitation  that  HONO  does, 
but  these  amounts  are  insignificant  compared  to  the  amount  of  energy  deposited  in  the  vibrational  modes 
(-50%). 
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Figure  4.  (a)  kTotal  and  fb)  kj  fsauaresl  and  kn  f circles]  and  (c)  branching  ratio  as  functions  of  enen 
Points  are  trajectory  results:  the  solid  lines  are  RRK  fits  to  the  trajectory  results. 
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Table  5.  Product  Energy  Partitioning* 


Reaction  I:  CH2NN02  ->  H2CH  +  ONO 


Initial  Energy15 

Available 
Product  Energy 

<^Trans> 

<^'Rot>'  H2CH 

<ERot> »  ONO 

62.1 

27.1 

4.4  (16) 

3.0  (11) 

■— 

73.6 

38.6 

6.9  (18) 

3.9  (10) 

WmMmmm 

85.1 

50.1 

9.0  (18) 

5.5  (11) 

4.6  (9) 

96.7 

61.7 

10.8  (18) 

7.1  (12) 

5.8  (9) 

108.2 

73.2 

12.0  (16) 

7.6  (10) 

119.7 

84.7 

15.7  (18) 

9.7  (12) 

7.4  (9) 

131.3 

96.3 

16.6  (17) 

9.9  (10) 

9.2  (10) 

Reaction  II:  CH2NN02  — »  HONO  +  HCN  (no  secondary  reactions) 


Initial  Energy15 

Available 
Product  Energy 

<ETrans> 

<ERot>,  HONO 

<Eroi>’  HCN 

62.1 

90.0 

15.2  (17) 

■HB 

8.3  (9) 

73.6 

101.5 

26.3  (26) 

12.9  (13) 

85.1 

113.0 

33.0  (29) 

15.9  (14) 

96.7 

124.6 

36.2  (29) 

7.1  (6) 

17.8  (14) 

108.2 

136.1 

39.4  (29) 

7.6  (6) 

19.4  (14) 

j  119.7 

147.6 

41.0  (28) 

8.8  (6) 

21.0  (14) 

131.3 

159.2 

46.1  (29) 

9.7  (6) 

22.6  (14) 

Reaction  III: 


CH2NN02 


HONO  +  HCN  — »  HO  +  NO  +  HCN 


Initial 

Energy 

Available 

Product 

Energy 

<^Trans>’ 

(HO+NO) 

+  HCN 

HCN 

<ETrans>-d 

HO+NO 

<ERot>*d 

HO 

<Erm." 

NO 

62.1 

38.0 

■ 

_ 

_ 

73.6 

49.5 

— 

— 

— 

— 

— 

85.1 

61.0 

5.3  (5) 

11.8  (10) 

19.4  (32) 

1.5  (2) 

96.7 

72.6 

9.2  (7) 

2.8  (2) 

26.1  (36) 

3.2  (4) 

7.6  (10) 

108.2 

84.1 

12.0  (9) 

7.1  (5) 

29.1  (35) 

6.2  (7) 

5.8  (7) 

119.7 

95.6 

11.8  (8) 

11.1  (8) 

27.0  (28) 

5.5  (6) 

11.1  (12) 

131.3 

107.2 

12.5  (8) 

13.1  (8) 

30.4  (28) 

5.5  (5) 

7.8  (7.3) 

a  Energy  given  in  kcal/mol. 

15  Includes  the  zero-point  energy. 

c  Percentage  of  available  product  energy  of  unreacted  HONO  +  HCN. 
**  Percentage  of  available  product  energy  after  HONO  decomposition. 


Average  product  energies  for  trajectories  resulting  in  secondary  decomposition  of  HONO  are  compared 
to  those  that  do  not  undergo  HONO  decomposition  in  Table  5.  The  percentages  of  available  product 
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energy  in  the  columns  labeled  <Errans>,  (HO+NO)  +  HCN  and  <%<„>,  HCN  in  Table  5  are  percentages 
of  the  available  product  energy  upon  HONO  formation  from  CH2NN02.  These  results  are  directly 
comparable  to  the  average  relative  translational  energy  and  average  HCN  rotational  energies  for  the  HONO 
trajectories  that  do  not  react.  The  rotational  energy  distributions  of  HCN  do  not  differ  significantly 
between  the  two  subsets  of  trajectories.  For  the  trajectories  that  result  in  secondary  decomposition  of 
HONO,  the  average  HCN  rotational  energy  is  slightly  less  than  that  corresponding  to  unreacted  HONO. 
The  average  relative  translational  energy  distributions  for  trajectories  resulting  in  secondary  decomposition 
of  HONO,  however,  differ  substantially  from  those  corresponding  to  unreacted  HONO. 

Figure  5  shows  the  relative  translational  energy  distributions  upon  HONO  formation  for  trajectories 
that  do  not  undergo  secondary  decomposition  (black  filled  bars)  and  those  that  do  (grey-shaded  bars)  for 
131  kcal/mol  total  energy.  The  distribution  of  the  trajectories  that  result  in  secondary  decomposition  of 
HONO  is  peaked  near  10  kcal/mol  and  has  an  average  value  corresponding  to  8%  of  the  available  product 
energy.  We  estimate  that  as  much  as  75%  of  the  energy  available  to  products  goes  into  the  vibrational 
modes  of  the  HONO  that  will  subsequently  decompose.  The  distribution  corresponding  to  unreacted 
HONO  peaks  at  55  kcal/mol,  and  has  an  average  value  of  29%  of  the  available  product  energy.  We 
estimate  that  only  50%  of  the  available  product  energy  goes  into  vibration  for  this  subset  of  trajectories 
resulting  in  HONO  formation.  It  is  clear  that  the  destiny  of  the  nascent  HONO  is  dependent  upon  how 
the  energy  released  in  the  primary  decomposition  of  CH2NN02  is  distributed  among  the  internal  and 
translational  modes  of  the  fragments.  Product  energy  averages  for  the  diatomics  resulting  from  HONO 
decomposition  indicate  that  both  HO  and  NO  have  little  rotational  excitation  but  substantial  translational 
excitation. 

Product  energy  distributions  for  reaction  (I)  are  unremarkable.  The  translational  energy  distributions 
are  peaked  at  small,  nonzero  energies  (<9  kcal/mol)  over  the  entire  energy  range  and  are  similar  to  those 
calculated  for  simple  bond  scission  reactions  in  other  systems  (see,  for  example,  Sewell  and  Thompson 
1991).  Rotational  energy  distributions  are  also  peaked  at  zero,  indicating  little  rotational  excitation. 

It  is  difficult  to  apply  the  results  of  this  study  to  the  experimental  measurements  of  Zhao,  Hintsa,  and 
Lee  (1988),  due  to  the  uncertainty  in  the  initial  state  of  the  nascent  CH2NN02  upon  primary 
decomposition  of  RDX,  and  the  ambiguity  in  the  assignment  of  the  parentage  of  the  resulting  final 
products.  Based  on  these  calculations  and  what  is  known  about  the  PES,  it  seems  reasonable  to  expea 
reaction  (I)  to  play  some  role  in  the  decomposition.  However,  Zhao,  Hintsa,  and  Lee  did  not  see 
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igure  3.  Relative  translational  energy  distributions  for  HONO  +  HCN.  Distribution  with  black  bars 
correspond  to  HONO  that  do  not  undergo  secondary  decomposition  and  grey  bars  correspond 
to  HONO  that  undergo  N-0  scission.  Ticks  along  the  abscissa  are  spaced  5  kcal/mol  apart. 

experimental  evidence  for  this  reaction  (Zhao,  Hintsa,  and  Lee  1988).  Rather,  the  experiments  have 
provided  compelling  evidence  for  the  concerted  reactions.  The  results  of  this  study  show  that  if  sufficient 
energy  is  partitioned  statistically  among  the  internal  modes  of  MN,  then  both  reactions  (I)  and  (II)  will 
occur,  with  reaction  (I)  being  more  probable.  Additionally,  the  decay  curves,  rate  coefficients,  and 
branching  ratios  show  no  anomalous  behavior.  One  explanation  that  would  satisfactorily  explain  the 
experimental  result,  without  contradicting  the  simulation  results,  is  that  the  internal  energy  of  the  nascent 
MN  formed  from  primary  decomposition  of  RDX  is  not  distributed  statistically  in  the  molecule,  but  placed 
preferentially  in  vibrational  modes  that  couple  strongly  into  the  concerted  molecular  elimination  routes. 
We  will  pursue  this  possibility  in  a  separate  study. 

5.  SUMMARY 

We  have  presented  molecular  dynamics  simulations  of  the  unimolecular  decomposition  of  CH2NN02. 
The  potential  energy  function  used  in  this  study  describes  two  competing  reaction  paths: 
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ch2nno2  ->  h2cn  +  no2 

and 

CH2NN02  ->  HCN  +  HONO.  (II) 

The  bond  dissociation  energy  for  N-N  bond  scission  in  (I)  is  35  kcal/mol;  this  reaction  has  no  back 
reaction  barrier.  The  activation  energy  of  (II)  has  been  estimated  to  be  31  ±  4  kcal/mol  based  on  the 
ab  initio  calculations  (and  corrections)  by  Mowrey  et  al.  1990.  The  transition  state  for  this  reaction  is  a 
five-centered  cyclic  structure  with  a  classical  barrier  of  37  kcal/mol.  The  thermal  activation  energy  for 
(II)  is  31.8  kcal/mol. 

Microcanonical  rate  coefficients  were  calculated  by  using  classical  trajectories  at  energies  between 
61  and  131  kcal/mol  (including  the  zero-point  energy,  27.5  kcal/mol).  The  results  show  that  the  dominant 
reaction  path  is  (I),  but  at  high  energies  (II)  is  competitive.  The  decay  curves  and  corresponding 
branching  ratios  of  the  decomposition  are  linear  and  are  well-described  by  the  simple  statistical  RRK 
model  (Robinson  and  Holbrook  1972).  There  does  not  appear  to  be  an  anomalous  energy  dependence  of 
the  rates  of  (I)  and  (II),  as  suggested  in  a  previous  study  (Mowrey  et  al.  1990).  HONO  formed  from  (II) 
can  undergo  secondary  reactions,  including  N-0  bond  scission,  H-0  bond  scission,  and  hydrogen  atom 
migration.  The  N-0  scission  and  hydrogen  atom  migration  appear  to  be  energy  dependent.  The  H-O 
scissions  are  infrequent  events. 

Product  energy  distributions  for  (I)  are  unremarkable  and  suggest  that  most  energy  available  to  the 
products  remains  in  vibration.  Product  energy  distributions  for  (II)  are  dependent  on  the  fate  of  the 
HONO  molecule.  For  the  subset  of  HONO  molecules  that  do  not  undergo  decomposition,  approximately 
1/3  of  the  excess  energy  is  deposited  in  relative  translation  between  HONO  and  HCN,  and  approximately 
50%  in  the  vibrational  modes.  For  the  subset  of  HONO  molecules  that  decompose  to  HO  +  NO,  less  than 
10%  of  the  excess  energy  is  deposited  in  relative  translation  between  HONO  and  HCN.  Because  the  HCN 
rotational  energy  distributions  do  not  differ  between  the  two  subsets,  we  conclude  that  HONO  molecules 
that  subsequently  decompose  have  approximately  25%  more  energy  partitioned  into  the  vibrational  modes 
upon  CH2NN02  decomposition  than  the  HONO  molecules  that  do  not  decompose. 

Possible  discrepancies  between  the  simulations  (which  indicate  that  [I]  is  the  major  decomposition 
pathway)  and  experiment  (which  find  concerted  molecular  elimination  pathways  to  dominate)  could  be 
due  to  an  anomalous  energy  distribution  in  the  nascent  CH2NN02  upon  RDX  decomposition.  We  will 
investigate  this  possibility  in  another  study. 
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